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POISSON FORMULA FOR RESONANCES IN EVEN DIMENSIONS* 

MACIEJ ZWORSKlt 

1. Introduction. We consider scattering by an abstract compactly supported 
perturbation in R™. To include the traditional cases of potential, obstacle and metric 
scattering without going into their particular nature we adopt the "black box" for- 
malism developed jointly with Sjostrand p3[ . It is quite likely that one could extend 
the results presented here to the case of non-compactly supported perturbation as 
well - see |2l| for a natural generalization of "black box" perturbations. 

We review now the basic assumptions. We work with a complex Hilbert space 
with an orthogonal decomposition 

(1.1) n = H Ro ®L 2 (M. n \B(Q,R )), 
and with an operator 

P : H — > H , self-adjoint with a domain T> C TL 

(1.2) l R "\ B (o,iio)£> = # 2 ( R " \ B(0, Ro)) 

1r™\B(0M o )P = - A|R"\B(0,i? ) , 

which satisfies 

(1.3) 3 k such that 1b(o,h ) (P + i) ~ k is of trace class , 

(1.4) P>-C, C>0. 

These assumptions guarantee that the resolvent R(X) = (P — A 2 ) -1 continues mero- 
morphically as an operator 7Y C omp — * l-'ioc from Im A < 0, A 2 ^ <t pp (P), to C when n is 
odd and to A, the logarithmic plane, when n is even. The poles of this meromorphic 
continuation are called resonances. At A ^ all reasonable definitions of multiplicity 
agree. We can for instance say that the multiplicity of a pole at A ^ is given by the 
rank of the polar part of 12(C) near A - see The situation is more subtle at and 
rather than go into a detailed discussion we will take the multiplicity required by the 
trace formula - see [^8] , , |3^] for the discussion of the resolvent near 0. 

If U(t) is the wave group for the operator P and Uo(t) is the free wave group we 
consider the natural wave trace: 



u(t) = tr [/(*) ~ U (t) , 

which is an even distribution in t 6 R. The notation used here is somewhat informal 
since U and Uq act on different spaces - see . The correct definition is given by 

(1.5) u(t) tr (U(t) - ^\B(0,R Q )Uo(t)i M "\B(QM )) +tr ^B(OM ) U o(t)^B(Q,R Q ) ■ 
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In odd dimensions the following Poisson formula was established in increasing degrees 
of generality by Bardos-Guillot-Ralston |Q, Melrose [^4[,p5| and Sjostrand-Zworski 
& 

(1.6) f 1+1 u(i) = f™ +1 ^m(A)e lA|t| , 

AeC 

ra(A) = multiplicity of A as a resonance of P , 

in the sense of distributions on R. The observation that we only need to multiply 
by t n+1 was made in |}5). The formula also holds exactly for super-exponentially 
decaying perturbations as was pointed out by Sa Barreto-Zworski [l9| . 
We note that for t > the trace formula is equivalent to 

(1.7) ^(A)=X>(C)^(A-C), 0eC c °°((O, TO )). 



The original proofs of (1.6) were based on Lax-Phillips theory |L2j and in par- 
ticular on the strong Huyghens principle. The extension to the case of hyperbolic 
surfaces by Guillope-Zworski Q provided a proof which does not require the strong 
Huyghens principle and is also applicable in the euclidean case . It is based on the 
Birman-Krein formula and "global minimum modulus" estimates on the scattering 
determinant. That was followed by a local trace formula of Sjostrand (2l) the proof 
of which did not involve any scattering theory but also used some "local minimum 
modulus" estimates for determinants of some holomorphic matrices. Sjostrand's for- 
mula specialized to the even dimensional compactly supported case gives the following 



weaker version of (1.7) 



(1.8) ^(A) = ^ m(C)<MA - C) + 0«A)-°°) , 

Q = [1/2, 3/2] + i[0, 1/2] , G C c °°((0, oo)) . 



We r emark however that the semi-classical local formula of 1 21 is much stronger than 
By using a "local minimum modulus" theorem in the argument of M, [p5| we 



can strengthen (L8) to obtain a global formula. This extension was motivated by a 



question asked by Vodev (see Sect. 3). 

Theorem 1. Let P be an operator sastisfying the assumptions (1.1)-(1.4) and 
let u(t) be its normalized wave trace given by ( fi-4 ). Let A p be an open conic neigh- 
bourhood of the real axis as shown in Fig.l, er(A) the scattering phase of P and let 
tp £ C C °°(R; [0, 1]) be equal to 1 near 0. Then 

(1.9) u(t) = ^ m(A)e 4A l t l + ^ m(A)e 



iX\t\ 



AeA p A 2 e CTpp (j=)n(-oo,o) 

Im A<0 



+ m(0) + 2/ ijj(X) — (X)cost\d\ + v p ^(t), t^O, 
Jo d\ 

v p4 , e C°°(M \ {0}) , d*v p ^ = 0(r N ) ,Vk,N, \t\ — ► oo . 

The scattering phase, ct(A), is a standard object in scattering theory - see (ijj] 
and reference given there for background information and for the discussion of the 
"black box" case. Here it is normalized so that the Birman-Krein formula holds - see 
below. 
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2. Proof of the trace formula. To prove Theorem 1 we identify the subset of 
A shown in Fig.l with C \ e" /2 S+, where e i7r / 2 R is the cut. Then the resonances are 
symmetric with respect to the cut and they coincide with the poles of the scattering 
determinant s(A) = det 5(A) . The unitarity of 5(A) for A > implies the usual rela- 
tion 5(A)- 1 = 5(A)* for Re A > 0. Hence s(A)" 1 = 1(A). To simplify the discussion^ 
we will consider the scattering matrix in A p n {Re A > 0} only and define s(A) in 
ApflRc A < so that the scattering phase defined by ct'(A) = (i /2tt)s' (X) / s(X) is even 



in A - see (|2_. 



The assumption ( |l.3| ) guarantees the existence of m such that for p > 

(2.1) m(X)<C p r m+e , Ve>0, 

|X|<r 

see [^5). This is deduced from the polynomial bounds of Vodev |30),|3l| which for 
p < ir/2 (all that is needed here) follow also from the earlier estimates of Sjostrand- 
Zworski ||. 
We now put 

(2.2) P p (X)= [] E(X/C,m) m ^ , E(z,p) = (l-z)e z+ ^ + - + ^ , 

CeA p \R 

where the bound ( [j.l| ) guarantees the convergence of the Weierstrass product. 

The poles of s(A) and s(A) _1 in A p coincide (with multiplicities) with the zeros 



x so that we do not need to consider global analytic properties of S(\) 
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of -Pp(A) and P p (— A) respectively. Hence we can write 

(2.3) »(A) = e B ' W ^J, Ae A p n{RcA>0}, 

where g p is holomorphic in A p n {Re A > 0}. We now extend g p (X), and consequently 
s(A), to A p by setting 

9 P {-X) = -9p(A) • 

That clearly implies that s'(X)/s(X) = s'(-X)/s(-X) for A e K \ {0} and further 
analysis shows that this identity holds through A = - see Sect. 3. 

We want to estimate the function g p . For that we need to estimate s(A) away 
from its poles and that is done exactly as in @],[^5| (see also Jl8| ). We write 

(2.4) S(X)=Id+A(X), 

A(X) = C n X n ~ 2 E^(-X)(I + K(X, Ao))" 1 ^, x]*^ 2 (A) , 

W : L 2 (R n ) — ► L a (S n-1 ) , W{0,x) = e iX < x > 9 '> p(x) , p e C~(R n ) , 

and where K(X, Xq) is the operator constructed in Sect. 3 of p3] ]: 

(2.5) K(X, A ) = [A, X o]i?o(A)(l - X i)x* ~ [A, X 2 ]i?(A )x 4 + X 2 (A 2 - A 2 )i?(A )x 4 , 

Xi G C C °°(]R") , xo = 1 near S(0, i? ) , Xi = 1 near supp Xi-i , 
R(X ) = (P-Xl)-\ ImAo «0, R (X) = (-A - A 2 )" 1 . 

To estimate s(A) we will first estimate ||(J + K(X, Ao ) ) 1 1 1 and that is based on the 
inequality 

(2.6) || ( / + /, ( A,A ))-||< det ( /+ '^^' m+1 ) 



det(/ + K(A,A )) m+1 )| 



from (j], Theorem 5.1, Chap.V. 

Exactly as in pj|, where we followed |p^|,p3]|, we see that for A 6 A p 

| det(J + K(X, A ) m+1 )| < Ce c|A| '" +5 . 

Using the lower modulus theorem of H. Cartan^J - see for instance Theorem 4, 
Sect. 11. 3 - we obtain a lower bound: 

\det(I + K(X,X a ) m+1 )\>Ce- Crm+ ^\ X € £>(r, pr/C) \(J D(Xj, rj ) , 5><r?r, 



uniformly as r — » oo. From this and ( |2 .6| ) we obtain, as in ]35[ , 

|a(A)|<Ce c ('- m+, /'») B , AGU^pr/COVy^.r,-), ^r^jjr. 

j i 

If we take i)Cl then for every r there exists r/2C < k(r) < r/C such that the circle 
A — r\ = k{r) does not intersect any of the excluded discs. Then using the standard 
estimates for Weierstrass products and the maximum principle we see that 



cxp 



g p {X)\\P p {-X)\ < Ce^ (m+ ' >n , A e A p , , < p' < p. 



2 We remark that a much cruder estimate would suffice here hut it is nice to quote the optimal 
result which is useful elsewhere in the theory of resonances - sec MM and Sect. 8. 
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We then conclude (as in the proof of Cartan's theorem or yet easier as in the proof 
of Hadamard's factorization theorem)^) that 

(2.7) \d k x g p (X)\<C\X\ im+ ^ n - k , AeV, < p" < pi , 

where the symbolic property followed from Cauchy's inequalities. 

We can now prove the Poisson formula. As in ||,fjj| the starting point is the 
Birman-Krein formula: 

(2.8) u{t) = a'{t) + ^ 2cos(iA) +m(0) , 

, - ■ A2 GCTpp (P)\{0} 

cr(A) = — logs(A) for A > , a(X) = -a(-A) for A < , 

where the Fourier transform is of course taken in the sense of distributions. For the 
"black box" perturbation the proof was given by Christiansen in Sect.l of but it is 
classical for all well known scattering problems. We note that our definition of s(A) 
implies that we can set 

Let tp be as in Theorem 1 and even. We can then write 
i s'{X) 



a'{X) = 



2tt s(A) 



(1 - iP(\))g' (\) + A (logP p (-A) - logP p (A)) 



^ (A) 7^y) +/ ^ (A) 



where f p ^ G C£°(R). 

The argument of || , |3j| now easily gives 

t 2 Pu(t)=t 2p ( e l ^m(X) 
\xeA p 

\ 

+ e im +m{Q) + ^a'{t) + {l^g p {t)+f p ^{t) 

A 2 e^ PP (-P)n(-oo,o) / 

ImKO / 

for p large enough. Since g p is a symbol on M, 

has the properties stated in Theorem 1 and this completes its proof. 

We remark here that a posteriori the bound on g p on the real axis has to be 



much better than the bound provided by the estimate (2.7): we know the strength 
of the singularity of u at t — and the bound on the number of resonances gives an 
estimate on the strength of the singularity of the exponential sum. Hence for elliptic 
perturbations where ra = nwe obtain 

\d k x g P (X)\ <C fe , e (l + |A|)"+ e - fc , Ve>0. 



3 We can use for instance Caratheodory's inequality - sec |27 
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3. Review of applications. The basic application of the trace formula is in 
obtaining lower bounds on the number of resonances from the singularities of the 
wave trace. The basic Tauberian lemma was given in Sjostrand-Zworski p4| and it 
was applied there to problems in odd dimensions. One of the applications of the local 
trace formula of Sjostrand ||l| was the extension of those results to even dimension - 
see Theorem 10.1 there. That becomes even clearer when we use the global formula 



(1.9). One of the interesting consequencies is based on the trace formula of Guillemin- 
Melrose [|: 

Theorem 2. Let P be the Dirichlet or Neumann Laplacian on a connected 
exterior domain R™ \ O where O has a smooth boundary. Suppose that the there exists 
a non- degenerate closed transversally reflected trajectory of the broken geodesic flow of 
M n \ O such that no essentially different closed trajectory has the same period. Then 
for any e > there exists C e > such that 

^{ro(A) : |A|<r + C,, | Im A| < e log |A|} > r/C e . 

For more applications we refer to and |2l]] , Sect. 10. 

The new trace formula allows also an easy extension of some of the results of Ikawa 
on the distribution of resonances for several convex obstacles to even dimensions. That 
is particularly interesting in dimension two where most of the numerical studies were 
conducted - see for instance || for the discussion of symbolic dynamics. We remark 
that it is rather clear that the results of (2l[] would suffice for this purpose but the 
global formula makes the applications even more apparent. As an example we give 
the modification of the result of Q : 

Theorem 3. Let P be the Neumann Laplacian on l n \ O, n > 2, and let us 
assume that O = U£=i 3 wh ere Oj are mutually disjoint strictly convex obstacles 
with smooth boundaries satisfying the following condition: 

convex hull (& k U £>i) n O m = V k ^ m j= I . 

Then there exists a > for which 

Y^{m(A) : | Im A| < a} = oo . 



The next theorem answers a question asked by Vodev |32| and does not seem to 
follow from the local trace formula: 

Theorem 4. Let P satisfy the assumptions (l.l)-(l.S) with n > 4 even. In 
addition let us assume that is not an eigenvalue or a zero-resonance of P. Then for 
any 7 > and k 



(3.1) 



u{t) - J2 ™(A)e lA l*l 

Im A<7 log I A| 



<C k r n+2 - k , t>t k > (n + k)/j. 



We remark that in odd dimensions the polynomial bound on the number of resonances 
and the global formula (L6) imply that the rig ht h and side of (3T) can be replaced 
by C(exp(— at)), a > 0. It is quite clear from (1.9) that to establish Theorem 4 we 
need to understand the behaviour of cr'(A) as A — * 0+. A finer analysis based for 
instance on [|l(J , |ll[] , should show that the estimate above is optimal and that in fact 
there exists an asymptotic expansion as t — > 00. Also, we did not attempt to study 
the more involved two dimensional case. 
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Proof. We start by observing that 
k 



fce 



J2 m(A)e lA|t| - J2 ™(A)e iA|t| I : C, 

. Im A<p|A| Im A<7 log |A| 

a k > , t > t k > (n + fc)/ 7 , 



which follows from the bound (2.1). In fact, the derivative of the difference can be 
estimated by 



-a'ut 



+ / x k X - lt dN(x) < e~< 1 + CM-7(M™+fc)/7) + Ct x --y(t-{n+k) /-,)-! dx 



t > (n + fc)/7 + e. 



(3.2) 



Thus we need to study the behaviour of 
S '(A) 



S (A) 



tr S{X)*S'{X) = tr (7 + A(X)*)A'(X) , 



as A — > 0+, where we used the notation of ( |2.4| ). 

The assumption that P has no resonance at zero implies that the cut-off resolvent, 
xR(X)x, X G C c °°(R n ), x = 1 near B(0,R Q ), is holomorphic in (A, A"" 2 log A) for 
A| < e. To see that we recall that the free resolvent, i?o(A), is of the form R' (X) + 
A"" 2 log AM(A) where R' (X) and M (A) are entire - see Sect.l of Q. Following [|||, 
Sect. 3, we write x-R(A)x = x(<3o(A)x + Qi(Ao)x)(/+-KXA, Ao)) _1 where Qo is a cut-off 
free resolvent and K is given by (2.5). Analytic Fredholm theory shows that when 
xR(X)x is bounded, it is a holomorphic function of A and A™~ 2 log A - see |2^], |3Cj 
for more details. 

To study (3.2) we need a different representation of A(X). If we recall the definition 
from Sect. 3 of |35|, A(X) comes from the radiation pattern of i?(A)(— [A, x]e^"'^). 
To obtain a formula similar to (2.4) but involving R(X) rather than (I + K(X, Ao)) _1 
we take Xj as in (p.5[) and write 



(1 - x 2 )i*(A)xi = R (X)(-A - A 2 )(l - X2)R(X) X i = Jfc(A)(-[A, xa]-R(A)xi) , 

since (1 - x 2 )(-A - A 2 ) = (1 - x 2 )(-P - A 2 ) and (1 - X2)xi = 0. This shows that 

A(X) = c n X n - 2 E^ (-A)[A, X2]R(X)[A, x]*E^ (A) , 

and the assumption on R(X ) im plies that A 2 ~" A(X) is holomorphic in A and A"~ 2 log A 
for |A| < e. From this and (^2) it follows that 

(3.3) ct'(A) = A™- 3 /(A,A"- 2 logA), A>0, 

with / smooth near 0. We then easily check that 

A"~ 3 .f (A, A"" 2 log X)tp(X) cos tXdX 
= (-l)*(n - 3)1/(0, 0) <-"+ 2 + 0(t- n+1 ) , t -> oo , 

which completes the proof as the estimates for derivatives clearly hold as well. □ 

Finally we compare this result with the estimates on the heat trace. As a conse- 
quence of well known estimates on heat kernels, Sa Barreto-Zworski p0[ showed that 
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when P = — A + V and P has no resonances with ImA < (that is no eigenvalues 
and no zero resonance) 

tr(e- tp -e tA ) =0(t~9 +1 ), t>0. 

Werner Miiller pointed out to the author that for the behaviour as t — ► oo it is more 
natural to study cr'(A) near A = using the heat version of the Birman-Krcin formula: 



POO 

(3.4) tr (e- tp ~ e tA ) = / e~ tA V(A)dA + V 



-tp,j 



where to make sense of the trace we used the convention employed in the definition 
of u, (1.5). Hence under the assumptions of Theorem 4 but for any n > 3 we obtain 
from its proof 

(3.5)tr(e- tp -e* A )= ]T e""* + ±T g - l) /(0, 0) + 0(t~t+4) , 

Mi <o 



where / is as in (3.3). In odd dimension it is a function of one variable, A, only. 
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